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ABSTRACT 

We briefly review a hamiltonian path integral formalism developed earlier by 
one of us. An important feature of this formalism is that the path integral 
quantization in arbitrary co-ordinates is set up making use of only classical 
hamiltonian without addition of adhoc Ti^ terms. In this paper we use this 
hamiltonian formalism and show how exact path integration may be done 
for several potentials. 
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1. Introduction 



The Feynman path integral formahsm has been in use for over forty years 
and has been successfully applied to a variety of problems^. An important 
development in the past fifteen years has been exact treatment of large num- 
ber of potential problems within the path integral formalism. Duru and 
Kleinert gave an exact path integral treatment of H- atom Green function 
using Kustaanheimo-Stiefel transformation^. The arguments originally used 
by Duru and Kleinert were formal and the manipulations used lacked math- 
emaical justification. The H atom problem was given a correct treatment by 
Ho and Inomata'* soon after Duru and Kleinerts work. It has since been pos- 
sible to give an exact path integral solution of several problems in quantum 
mechanics using the new technique scaling of local time. Notable among these 
are Morse oscillator^, Rosen-Morse oscillator^, Poschl-Tcllcr potential''', Hart- 
mann potential^, Hulthen potential*^, 6 - function potential^*', square well^^, 
and many other potential problems problem^^'^^. Besides the scaling of local 
time in combination with change of variables, the techniques of adding new 
degrees of freedom has been useful in most of the above mentioned applica- 
tions. 

The path integral formalism can be derived from quantum mechanics and 
this aspect has been a subject of extensive study in the hterature. One can 
arrive at different types of path integral representations of the propagator by 
working with different bases. Thus, for example, we have configuration space 
path integrals, phase space path integral or coherent state representation of 
path integrals. In all these derivations quantum mechanics is assumed. 

One of the most important features of the path integral formalism is that 
it gives an alternative route to quantization. In use of the path integrals as 
a scheme of quantization a prescription for setting up path integral is as- 
sumed and quantum mechanics is derived. The path integral quantization 
can be set up starting from classical lagrangian or hamiltonian form of the 
classical action for the system of interest. In the hamiltonian formalism of 
classical mechanics one can write dynamical equations in any pair of conju- 
gate variables related to a given set by a canonical transformation. However, 
it appears very difficult to implement general canonical transformations in 
quantum mechanics even though there are a large number of investigations 
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on this subjcct^"^. However, the point transformations do not present any 
difficulty in quantum mechanics. So it is of interest to ask if path inte- 
gral quantization scheme, based on lagrangian or hamiltonian form for the 
action, can be formulated in arbitrary co-ordinates. The answer has been 
known in the lagrangian formulation known for a long time. The configura- 
tion space path integral can be set up in such a way that it works in arbitrary 
co-ordinates utilizing the configuration space form of classical action^^. 

In a quantization scheme developed by one of us^^'^^, quantization is car- 
ried out in arbitrary co-ordinates using the hamiltonian form of path inte- 
grals. A unique feature of this scheme is that addition of 0{h^) terms is not 
necessary in any co-ordinate system and only the classical hamiltonian is used 
to set up hamiltonian path integral quantization in arbitrary co-ordinates. 
This formalism makes an essential use of the idea of local scaling of time to 
introduce a path integral representation for quantum mechanical propagator. 

In ref. 17 the basic property of the canonical path integrals, that only 
classical hamiltonian is needed for quantization in arbitrary co-ordinates, was 
established; discussion of further properties and also of applications was not 
taken up. In ref. 18 some important properties of the canonical path integral, 
such as scaling of local time, equivalence with lagrangian form of path inte- 
gration have been investigated in detail. In this paper possible applications of 
the scheme were indicated only briefly. In the past fifteen years many exactly 
solvable quantum mechanical problems have been treated by means of path 
integration by relating the given problem to another problem for which ex- 
act path integration can be done or is already known. In this paper we show 
that our scheme offers results which can be used in a very simple manner to 
relate different path integrals. Although we, generally, follow the techniques 
already available in the hterature to give exact treatment, the details of the 
treatments are different. There are some new points which require special 
and careful treatment. For example, initial condition for the propagator has 
to be treated carefully. 

In the next section we summarize all the necessary results on hamiltonian 

path integration from ref. 18. Some aspects of changes of variables and the 
technique of adding new degrees of freedom are discussed in Sec. 3 and Sec. 4 
respectively. Application of the hamiltonian path integral method to exact 
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solution of the quantum mechanical problems is discussed in the Sec. 5 to 7. 

2. Definition and properties of the hamiltonian path integrals 

In this section we briefly recall the definition of the hamiltonian path 
integral^^'^^. The proofs and detailed discussion of the results can be found 
in our earlier paper^^. We shall assume that all the hamiltonian like func- 
tions appearing in this paper are independent of time and are quadratic in 
momenta. 

2.1 Definitions : Given a phase space function H{q,p), we introduce 
a quantity {qt\\qoto) which stands for the short time approximation and in 
terms of which a path integral for finite times is constructed. We assume 
(O'^lko^o) to be of the form 



The definition given below is such that (0) satisfies the semi- group property 
with respect to the measure p{q)d^q 



when the integral is calculated in the stationary phase approximation. The 
'mixed short time propagators {qt\piti) and (j9iti|go^o) are defined below. Let 
7i and 72 be two classical trajectories with boundary conditions as indicated 
below. 




(1) 





7i : T 



{q{T),p{T)), for to<r<t 



1 



(3) 



^(^o) = go, Pih) = Pi 



(4) 



72 : T 



(g(r),p(r)), for h<T<t 



(5) 



P{ti)=Pu qit) = q 



(6) 
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We define 



{qt\p,h) = {2nh)'-/'JD^exp[zS++{qt,p,h)/h] (7) 



{Pitilqoto) = (27r/i)'"/V^~ exp[«5__(piti, qoto)/h] (8) 



where 



and 5"++, 5" are Legendre transforms of the classical action along the two 

trajectories, 

S++{qt,Pih)=puq[ + I {p,dq {r)-H(q{r),p{r)))dT (11) 
S-4qoto,pih) = -puq[ + / {Pidq' - H {q^r) , p{t))) dr (12) 

where 



qi = qih), qi=q (h) (13) 

With the above notation we now introduce the definition of the short time 
propagator (STP) {qt\\qoto) by 

{qt\\qoto) = I ^ / rf> {qt\piti) ipih\qoto) (14) 
VP(?)p(?o) 



In (^ to (|T^ it is understood that the independent variables are those 



explicitly shown on the left hand side. It is also understood that that the 
right hand sides have been expressed in terms of the independent variables 
using the equations (||) , (H) for the classical trajectories 71 and 72. The 
function 5*++ is the generator of canonical transformation connecting 'old' 
co-ordinates and momenta at time t to the 'new' ones at time ti. Similarly, 
S is the generator of canonical transformation connecting the conjugate 
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variables for the trajectory 71 at time ti to those at time to- It must be 
emphasisized that S±± are need for short times only. 



We shall define two path integrals, the first one is without any scaling 
of time and the second one with scaling of time. The first hamiltonian path 
integral K[H, p] defined as a summation over histories with (g'e||goO) inserted 
for short times: 



with e = t/N and = q. 

Motivated by the importance of the local scaling of time in the path 
integral formalism, we now introduce the second hamiltonian path integral 
which will be described as hamiltonian path integral with scaling of time. The 
definition makes use of hamiltonian path integral HPIl. A special case of the 
hamiltonian path integral with scaling was found suitable for quantization in 
arbitrary coordinates^'^. 

Given a hamiltonian function H{q,p), measure p{q), and a strictly posi- 
tive function a{q), we introduce second hamiltonian path integral with local 
scaling of time, to be denoted by IC[H, p, a] as 

}C[H,p,a]{qt;qoO) 



a{q){H — E) as 'hamiltonian' function. T/ie two hamiltonian path integrals K 
and K. will henceforth be called HPIl and HPI2 respectively . This completes 
the definition of HPI2. We shall now give several remarks concerning the 
definintion of HPI2 which will be helpful in clarifying our approach. 

1. Our first remark concerns the definition of HPI2 and use of local scaling 
of time in this paper. Although the definition of HPI2 is inspired by the 
expression and formulas appearing in the literature, we do not utilize 
or rely on any of the other details, mathematical or otherwise, from 
papers in the existing literature on local scaling of time. 



K[H,p]{qt;qQto) = lim / [| p{qk)dqk / W {qk+i4<lk^) (15) 
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2. The second remark is about realtionship between HPIl and HPI2. The 
path integral HPI2 generahzes HPIl and conversely HPIl becomes a 
special case of HPI2 in the sense that for trivial scaling, when the 
scaling function a{q) is a constant independent of time, HPI2 becomes 
equal to HPIl. For this reason we shall continue to describe HPI2 as 
the hamiltonian path integral with scaling. 

3. The third remark is on the appearance of the functions a{q) and p{q) in 
the definition of HPI2.As far as the definition is concerned the HPI2 is 
defined for arbitrary scaling function a{q) and measure p{q) just as the 
hamiltonian function H{q,p) is kept arbitrary in the definition of any 
path integral. Specific values such as those given by a{q) = y/g and 
p{l) = y/9 give rise to the path integral scheme suitable for quantization 
in arbitrary coordinates ( see Sec. 3.2 below). 

4. Our last reamark concerns the properties and applications of HPI2. 
The results on scaling developed earlier^^ and others obtained here have 
been applied to exact path integration by relating the path integral 
for the problem to be solved to another problem for which answers 
are already known or can be worked out. When one attempts to do 
this one has to, in general, deal with other intermediate expressions 
involving HPIl or HPI2 with different choices of a{q) and p{q). It is 
for this purpose a{q) and p{q) have to be kept arbitrary while defining 



In the remaining part of this section we shall summarize important properties 
of the two hamiltonian path integrals 

2.2 Descretized form : In the limit e — 0, the only important terms are 
are the terms of order e and only these need be retained in the definition of 
STP. Retaining such terms only, the canonical STP is equivalent to 



HPI2. 




exp \ ipi{q' - go)/^J ^ 

r _ e d'^H{q,p) ^ e d^H{qo,p) 
2 dq^dpi 2 dq^dpi 



(17) 



Inserting (|17D in ([15|) 
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N-1 \ N-1 



K[H,p]{qt;qoto) = lim / [| p{qk)dqk] Y[{qm+i4QmO) (18) 

\k=l I m=0 

gives discrete form for the first hamiltonian path integral. Notice that prod- 



ucts over k and over m appearing in ([ISD have different ranges. Hence this 



expression involves — 1 fold q integrations and A^— fold p integrations. It 
is easy to do the p— integrations in the STP if the hamiltonian is quadratic 
in momenta. This gives a lagrangian form of STP which can be brought to 
a standard form making use of the McLaughlin Schulman trick. This makes 
a connection of our hamiltonian path integrals with lagrangian form of path 
integral possible 

2.3 Immediate consequences of definition of path integrals : We make 
some useful observations which follow directly from the definitions and the 
discrete form for the short time approximation. 

(a) Addition of a constant to hamiltonian H is equivalent to multiplying 
the HP II by a phase 

K[H - E, p\{qt;qoO) = exp{tEt/h)K[H, p]iqt;qoO). (19) 

(b) The dependence of the HPIl and HPI2 on the measure p is particularly 
simple. Thus for example we have 



p{q)p{qo)K[H,p]=K[H,l] (20) 



K[H,p,] = ^{a{q)a{qo))K[H,p2] (21) 
where pi and p2 are related by 

P2(g) = a(g)Pi(g). (22) 

Similar relations are true for HPI2 JC also. 

(c) When the scaling function a is a constant, say c, we have 

lC[H,p,a = c] = K[H,p] (23) 
This equation is equivalent to 
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/dE f°° 1 

j^expi-iEt/h) daK[ciH - E),p]{qt;qoO) = -K[H,p]iqt;qoO) 

(24) 

(d) In the special case when the function H{q,p) is independent of q^ for 
some j and depends on pj alone, the HPIl depends only on the difference 

(g — go), apart from an overall factor of \fp{^)p{(lo) ■ This is easily seen by 
doing the corresponding p integrations after writing the discrete form for 
HPIl. 

2.4 Schrodinger equation for /C : If we take 

H = ^g''PiPj + V{q) (25) 

the HPI2 IC[H, a, p]{qt; go^o) as defined in (jl^) satisfies the Schrodinger equa- 
tion 

d 

ih—}C[H, a, p] = {Ha)op}C[H, a, p] (26) 

with 

(H.u = -^p-'"<^r'^" (l^^^'^^ly) p"''"^^^r'^" + n^) (27) 

and /C has initial value given by 

hm 1C[H, ^, ^]{qt; goto) = (l/p(g))5(g - go) (28) 

For the special case when a = p = ^ the HPI2 IC[H, y/g]{qt; goto) 
with scaling, as defined in (0) satisfies the Schrodinger equation 

^fl^/C[i^, v^, = HIC[H, ^, ^ (29) 

with 
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and K, has initial value given by 



lim 1C[H, ^]{qt; goto) = 9-'^\q)S{q - go) (31) 

The relations (|26|) to (0) are the central results in our scheme. In our 
scheme JC[H, ^/g, ^/g] is the candidate for the quantum mechanical propaga- 
tor for a Hamiltonian of type (p5[) because the right hand side of (pOD can be 
recognized as the Schrodinger operator —(h /2m)V^ + V{q). 

The proof of ( P7| ) and ( PB[ ) is long and can be given in several ways. Here 
we indicate, in brief, one of the methods used in ref. 18. For this purpose 
consider the "propagation of functions of q in time" governed by the HPI2 
K. and define 

^{q,t) = Jp{qo)dqo}C[H,p,a]{qt-qoO)^{qo) (32) 

The quantity \E'(g, t + At) — \E'(g, t) is computed up to the first order terms for 
small At. To do this we by first insert discrete form for the HPIl K[a{H — 
E),p\ in the right hand side of (|16D and use the resulting expression for 



1C[H, p,a{q)] in the right hand side of (|3^) . Because in (|T6]) integral of 
K[a{H — E), p]{qa,qQO) over a is needed, it is not possible to replace K 
by a single STP and, we have to retain full discrete expression for finite A^. 
The expression (|32D involves fold g— and fold p— integrations as well 
as integrations over E and a. The integration over E results in a Dirac 
delta function implying a linear constraint in a and At with functions of g 
appearing as coefficients. The a— integration is easily done using the delta 
function. The resulting expression is expanded in powers of At retaining all 
0{At) terms. The result of the remaining fold g— and p— integrations in 
( |32D can be expressed as an operator acting on \E' after long algebraic manip- 
ulations, he result can also be proved by making use of the correspondance 
of the HPIl with lagrangian form of path integral mentioned earlier at the 
end of Sec. 2.2. 

As a by product of the above result we also get the scaling formula given 
below. 

2.5 Scaling formula : The HPI2 IC[H, p, a] with scaling function a{q) and 
density p, as defined by ([161) above, is related to the HPI2 with trivial scaling 
and hence to HPIl with hamiltonian H — U{a) 

IC[H, p, a] = IC[H - U{a), p, 1] = K[H - U{a),p] (33) 
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The precise form of U{a) depends on H used. For H{q,p) given by (|25| ) we 
have 

U{a) = -— (g'^dSna)d,{lna) + 2d,{g'^d,{lna))) (34) 

Redefining H in the scahng formula (0), it can be put into a useful form as 
a relation between two HPIl. 

K[H,p]{qt;qoO) 
= IC[H + U{a),p,a]{qt-qoO) 

I t dE t°° 

= ya{q)a{qo) J ^^^^^ exp{-iEt/h) daK[a{H - E + U{a)), p]{qa; qoO) 

■ ■ ■ (35) 

This formula expresses HPIl without scaling in terms of a HPIl with scal- 
ing. The factor \J ot{q)oi{qo) can be absorbed into change of measure from p 
to p/a for the HPIl in the right hand side of (^) 



K[H,p]{qt-qoQ) 

/dE r°° 
j^exp{-iEt/h)J^ daK[a{H -E + U{a)),{p/a)]{qa;qoO) (36) 

Other useful forms of the above relations are obtained by integrating over 
time or taking inverse Fourier transforms. As an example, we can obtain 

J dtK[a{H-E),p]{qt-qoO) 

= , / dtK[{H-E + U{a)),p]{qt;qoO) (37) 

^Ja{q)a{qo) •'^ 

by integrating the scaling relation (|33[) over time. 



3. Point transformations and canonical path integration 

3. 1 Initial condition for canonical path integration : For a path integral to 
give rise to correct propagator, it is important that not only the correct 
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Schrodinger equation be reproduced, it must also satisfy boundary conditions 
appropriate to the co-ordinates chosen. In this section we will at first take 
an illustrative example and then discuss the general case of initial condition 
for the canonical path integration. As an example of initial condition on the 
propagator in arbitrary coordinates, consider the change from the cartesian 
co-ordinates to polar co-ordinates in a plane. 

X = r cos 6, y = r sin 6 (38) 

The right hand side of (|38|) is invariant, separately, under the following trans- 
formations. 



(i) r ^ r, 6^9 + 2m7r; (39) 

(ii) r ^ -r, 9 ^ 6 + {2m + 1) (40) 

Thus we have the relation 



-| oo 

5{x-xo)5iy-yo) = - E {5ir-ro)5{9-9o+27im)+5ir+ro)5[9-9o+i2m+l)7r]) 



m=—oo 



(41) 

Quantum mechanical propagator, {ft\fQto), in two dimensions satisfies 
the initial condition in cartesian co-ordinates 



lim(r,t|foto) = 5{f- fo) (42) 

Therefore, in plane polar co-ordinates the propagator, at equal times, reduces 
to the right hand side of (|4lD . Conversely for the change of variable from r, 9 

to X, y 

r = (x^ + ^/2)^/^ 9 = td.n-\y/x), (43) 

let us assume that a hamiltonian path integral, say HPIl, is constructed to 
satisfy the initial condition 

K[H, r] {r9to; ro^oto) = - ro)6{9 - 9^) (44) 

Because the right hand side of (|43|) does not change when (x, y) —>■ {—x, —y) 
, therefore, the equation 
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-5{r - ro)5{9 - Oq) = 5{x - Xo)5iy - yo) + + Xo)5{y + yo) (45) 
r 

correctly gives the boundary condition (^) in the cartesian co- ordinates. In 
general let g,p — >■ Q, P be a point transformation. 

Q = Q{q), P = Piq,p) (46) 
with the inverse transformations written as 

q = qiQ), p = p{Q,P) (47) 

Let h{q,p) and H{Q,P) = h{q{Q),p{Q, P)) be hamiltonian functions in the 
two sets of co-ordinates. Let us set up HPI2 in the two sets of co-ordinates q 
and Q directly using the scaled canonical path integrals using the same clas- 
sical Hamiltonian but expressed in terms of the appropriate canonical vari- 
ables. The two HPI2 /C[/i, y^^, v^J(gt; goO) and K,[H, ^^]{Qt-Q^Q) 
satisfy the same Schrodinger equation and by construction they obey the 
following boundary conditions in the limit t — 0. 

V9v VaM = 0; goO) = gl^\q)6{q - go) (48) 

IC[H, ^^]iQt = 0; QoO) = gl^\Q)SiQ - Qo) (49) 

Whenever the function q{Q) can be solved for Q uniquely, }C[H, 
will be obtained from }C[h, ^/g^] of (^) by expressing g in terms of 
Q. This relationship will not hold whenever the relation g = q{Q) cannot 
be inverted to give Q uniquely in terms of g. This is due to the fact that 
in such a case the right hand sides of (|4^) and (|49|) are not equal. Let 
A be a transformation on Q's such that q{AQ) = q{Q) . The set of all 
such transformations form a group. Following Moshinsky et al. we call it 
ambiguity group of the transformation q ^ Q, and denote it by A. Then 

S{q-qo) = \J\T.^iQ-^Qo) (50) 

In the above equation Qq stands for any one solution of the g = q{Q) and J 
is the Jacobian of the transformation q ^ Q . The relation between the two 
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propagators p8| ) and ( ^91 ) can the be written down 



= J2 ^[H. V9~2,V9^{Qt; AQoO) = ^ }C[H, y^, ^{AQt; QoO)(51) 



In the equation (|5lD it is understood that the variables q, go are to be ex- 
pressed in terms of Q and Qo by using the equations for the point transfor- 
mation. The resuh (^) follows from the fact that both sides obey the same 
Schrodinger equation and the same boundary condition. 

3.2 Quantization in arbitrary co-ordinates: We are now in a position to 
write the path integral for quantization in arbitrary co-ordinates needed for 
later applications. We are interested in the potential problems of quantum 
mechanics. The classical hamiltonian, hd, for these problems is of the form 

Ki=py2m + V{r) (52) 

where r denotes cartesian co-ordinates and p are the conjugate momenta. Let 
Q be any other set of co-ordinates related to r by a point transformation . 
Then the quantum mechanical propagator {ft\foO) satisfying the Schrodinger 
equation and the correct boundary condition 

lirn {rt\r qO) = 6{r — Vq) (53) 

can be written in terms of HP12 set up in Q, P variables as 
(ft|foO) 

A€A 

dE 



f fj fij too 

{9{Q)9mf" E / 7^exp(-zi?t/^) / daK[HE^^ma-Q,Q) 

/dE 
— — eM-^Et/n) / daK[HE, p = I] (AQa; QqO) (54) 



where He = ^/g{Hcl — E). The corresponding energy dependent Green 
function is defined by 

roo 

G{ft,fo\E)= dtexp{iE+t/h) (ftlfoO), (55) 
Jo 
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where Ej^ = E + iri,ri > and Mmrj ^ is taken at the end. For the energy 
dependent Green function we have the following path integral representation 
in arbitrary co-ordinates. 



G{f,ro\E) 

/•CO 

= {9{Q)9{Q^)Y"T. daK[^{H,i-E),^{AQaQoO) (56) 







J2 / daK[^{H,i -E),p= l](Aga,goO). (57) 



In the above Hd is the classical hamiltonian written in terms of Q and P . 
The expressions ( |5^ j and ( |5^ j define, respectively, the scheme for obtaining 



the quantum mechanical propagator and the energy dependent Green function. 
They involve the classical Hamiltonian only which is to he expressed in terms 
of appropriate co-ordinates and conjugate momenta, no ad hoc addition of 
0{h^) terms is necessary in any co-ordinate system. 

4. Relating path integrals 



In many of the exact evaluations of the path integrals in literature, the re- 
sults are be obtained by relating the desired proppagators to one for which 
answers are known. There are some notable investigations in which the au- 
thors actually compute the multi-dimensional integrations to complete path 
integration. Besides the use of local rescaling of time, the technique of adding 
new degrees of freedom with known, preferably trivial, dynamics has proved 
to be useful in relating path integral representations for quantum mechanical 
propagators of different quantum mechanical problems. Adding new degrees 
of freedom is like inverse of separating the variables in the differential equa- 
tion approach. We shall briefly recall this technique here. 

Let Hext{(l,P, P) be independent of Q and H{q,p) be such that for P = 
constant, say u , we have 

H,,t{q,p,P)\p=, = H{q,p) (58) 

We will say that Hext is an extension of H by addition of a new degree 
freedom Q, P with trivial dynamics. Let K and K'^^^ be the two HPIl using 
the unit measure and with H and Hext as hamiltonians. Then 
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K[H, l]iqt; qoO) = J d{/\Q) exp{iuAQ)K[He.t, l]iqQt; qoQoO) (59) 

Notice that Hext is assumed to be independent of Q. It, therefore, follows 
that K[Hext, 1] depends only on Q — Qo{= AQ) . The limits of integration 
over AQ are usually from — (X) to +00. However, if K[Hext, 1] is periodic in 
AQ the region of integration can be restricted to one period for the case 
u = integer. Similar results have been used in literature and one can easily 
arrived at the result by considering the differential equation satisfied by the 
the two path integrals appearing in (^). 

In this paper almost all the examples of path integrations to be discuused, 
with the exception of the hydrogen atom in the parabolic co-ordinates, the 
desired propapgator or the energy dependent Green function will be obtained 
by relating the corresponding path integral representation to another path 
integral which is known. A careful distinction has to be maintained between 
the path integrals which are related to quantum mechanical answers for prop- 
agator of some problem and other path integrals which appear only in the 
intermediate steps. For any quantum mechanical problem the propagator is 
given by (a) the corresponding HPIl in cartesian co-ordinates or an equiva- 
lent HPI2 expression with trivial scaling, or (b) expression (|5^ ) appropriate 



to the non-cartesian co-ordinate system selected, or (c) the expression in the 
line just before (|5^). 

Apart from the addition of new degrees of freedom, properties about of 
HPIl and HPI2 when p or a are changed will be used in obtaining relations 
between different path integrals. 

A few remrks about the notation to be followed are given in the following. 
We shall always use the bold face symbols Ko(gt|go^o) and Go{q, qo\E), with 
a zero subscript, to denote the full quantum mechanical solution for the prop- 
agator and the energy dependent Green function, respectively. Similarly, the 
hamiltonian for the problem to be solved will have a subscript 0. 

Different superscripts or subscripts on roman, script or bold face symbols 
will be used to denote other intermediate or related expressions . Many of 
the exact treatments in this paper begin with a series of observations and 
from these observations we "zero in" to one of the two required functions 

Ko(gt|go^o) or Go(g,gol^)- 

In the following sections the symbol ki(xt|xoO) will be used to denote 
harmonic oscillator propagator in one dimension for mass fi and frequency a; 
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kf'=(xt|xoO) = — exp (i ^ [(x^ + xl) cosut - 2xxo] ] (60) 

2m sm ut \ 2 sm ut J 

The propagators for harmonic oscillator in n dimension will be product of n 
such factors and will be denoted by k°'^'^ . Also in the following we shall use 
same symbol to denote the classical hamiltonian in different co-ordinates. 

5. Potential problems in one dimension 



5.1 The ax"^ + h/x"^ potential on half line: The solution for this problem 
will be obtained by relating it, following Duru^^, to two dimensional oscillator 
problem in a half plane. The derivation is written as a sequence of steps 
consisting of several small intermediate results and observations. 

Stepl : The hamiltonian for the problem given by 

2 

Ho = — + ax^ + hlx^. (61) 
2m 

The required propagator Ko(a;t|xoto) has initial value 

lim Ko(xt|xoto) = ^{.x — xq) (62) 

The propagator Kg , in terms of the hamiltonian path integral representation, 
is given by 

Ko(xt|xoO) = 1C[Hq, p = 1, a = l](xt; XQt) = K[Hq, p = l]{xt; xqO) (63) 



This is due to the fact that HPIl in (|63|), the last expression, is the right an- 
swer in the Cartesian co-ordinates. It obeys the correct Schrodinger equation 
and the correct initial condition at t = 0. 



Step 2 : The classical hamiltonian Hq looks like the hamiltonian H 



with po taken as suitable constant. Under a point transformation 
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xi = X cosh 9, X2 = X sinh 9 



(65) 



dxidx2 = xdxd9, ^/g = x (66) 

the hamiltonian Hi takes the form of harmonic oscillator hamiltonian in two 
dimensions 

for which the the exact quantum mechanical propagator is known. 

Step 3 : We set up the path integral quantization scheme for oscillator 
problem in the two sets of co-ordinate systems, viz., the plane polar co- 
ordinates (x, 6') and cartesian co-ordinates x = {xi,X2) and define 



/Ci = }C[Hi,l,l\{x,t;xoO) = K[Hi,l\{x,t;xQ) (68) 



and 



where 



IC2 = IC[Hi,^,^]{x,9t;xo9oO) (69) 

dE f 
—— expi-iEt/h) j daK{H2, ^{x, 9<7; a:„9„0) (70) 

H., = .(H.-E)=.(§-^jl-,^a.^-E) (71) 

The two functions /Ci and /C2 become equal when expressed in terms of either 
co-ordinates because the mapping xi, X2 to x, 9 is one to one. They are also 
equal to the quantum mechanical propagator in two dimensions for (p7|). 

/Ci = /C2 = kf ^(xit|xioO)kf ^(x2t|a;2oO) (72) 
k°^'^(xt|xoO) appearing in the above equation is the propagator for one di- 



mensional oscillator with mass m and frequency cu given by \j2a/m. 

Step 4 '■ Although Hi reduces to Hq for pg = const, but x dependences 
of the normalizations of /C2 and Kq do not match: 

Ko(a;, t = 0; xq, 0) = 6{x — xq) (73) 
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and 

/Cafx, e,t = 0; xo, Oq, 0) = -S(x - xo)5(e - ^o) (74) 

X 

The steps to relate these two have to take into account of this difference in 
the initial conditions at t = 0. In order to have correct initial value for the 
two dimensional problem, we define 

K{xt\xoto) = —^=Ko{xt\xoto) (75) 
y/xxo 

which obeys 

K{xt\xoto)\t=to = -5{x - xq) (76) 



Step 5 : Next We insert the HPIl of (|63D for Ko(a;t; XqO) in the right hand 
side of ( [75| ) and use the formula (p7|). Dividing (|35| ) by \Ja{q)a{qo) and using 
it we get 

K{xt\xoO) 

/roo 
dEexp{-iEt/h) da K[H3, p = l]{x, a; xqO) (77) 

where 

Hs = x{Ho + U{x)-E) 

1 2m ~'~ ( ~'~ 2m) ~'~ ^ } 

Step 6 : This problem of obtaining HPIl for in ( [76D is related to the 
HPIl for the extended hamiltonian H2 

^-K^^^^"'^"^) '''' 

for pg = {2mb/h +1/4) . Thus using (4.2) we get the relation 

K[H3,p= l](x,a;xoO) = J d{Ae) exp{iuAe)K[H2, p = v^](x0(t; xqM) 

(79) 
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where A6 = 6 — 60 and 



1-, 



2mb 1\ , , 

+ 7 (80) 



Stt'p 7 : From ( |7^ Kq is given in terms of Kq which can be expressed in 
terms of K\H2, 1] using (|77D and (^). Thus we obtain 

Ko(xt|xo,0) 
= ^xxqK. < xt\xoO > 

//•oo 
dE exp{—iEt/h) j daK[Hs, p = l]{x, a; xqO) 

d(Aa) exp^iuO) J dE exp^—iEt/h) J daK[H2, x] (x^a; Xq^qO) 

d{Aa) exp(zz/0)kr (81) 

Here the last step comes from equahty (|72|) aheady noted and is the 
propagator for the two dimensional problem [67|).We get the desired result by 
and following ref. 19 for the remaining steps and integrations. Omitting the 
details, the final answer is 



Ko(xt xoO) = — —I^ . ■ ^ exp 



f imtu \ 2 



(82) 



5.2 Morse Oscillator : We shall briefly indicate how the propagator for 
this problem is related to that for ax"^ + b/x"^ potential on half line. The 
Morse oscillator hamiltonian 

Hq = h y4exp(— 2ax) — B exp(—ax) (83) 

2m 

after a point transformation to a new coordinate r = exp{~ax/2) becomes 

Ho = —a\y, + {Ar^ - Br^) (84) 
2m 

Quantization of the Morse oscillator in the new variable r requires us to set 
up 
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I [ dE f°° 

= ^Ja{q)a{qo) J j^^exp{-iEt/n) d(TK[HE, p]{ra;roO) (85) 

where a = ^/g = (2/ar) and He is given by 

He = {2/ar){Ho-E) 

= (ar/2)(^ + -(Ar'-B) ^ 



2m o?r'^ ) 

= {ar/2)H2 (86) 

where H2 is the function in the brackets in the second hne of equation (|86|). 
The HPIl in r.h.s. of ( ^5]) can be related to HPIl for H2 using the scahng 
formula with a U{ar/2) appearing in the Hamiltonian. The HPIl for the 
resulting hamiltonian H2 — U{ar/2) is just the HPIl for the radial oscillator 
and is known from the results in the previous subsection. 

6. Coulomb problem in two dimensions 

The Coulomb problem in two dimensions is solved easily by our method in 
the parabolic co-ordinates as it gets related to harmonic oscillator in two 
dimensions. The steps leading to the solution are as follows. 

Step 1: The classical hamiltonian for the Coulomb problem in two di- 
mensions is 

Ho = ^-- (87) 

2m r 

P= iPx,Py),r = {x,y) (88) 
In the " 2-dimensional parabolic coordinates " ui,U2 defined 

X = ul — ul; y = 2uiU2 (89) 

the classical hamiltonian takes the form 

^0 = ^-4 (90) 

dxdy = Au '^duidu2, g^^^ = (91) 

Step 2 : For quantization of the Coulomb system in the cartesian, recall 
for cartesian co-ordinates the desired propagator Ko(rt|roO) is just the HPI2 
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with a trivial scaling, which is also equal to HPIl for the same hamiltonian. 
Thus we have 



Ko(rl|r-t)) = /C[i/o, 1, 1] = K[H,, 1] (92) 

Step 3: Quantization in parabolic coordinates ui, U2 will proceed via path 
integral HPI2 with hamiltonian Hq, scaling function a = Au^. Thus we should 
set up /C[i/o, 4m^, 4u^]('ut; uqO) 

/dE r°° 
——exp{-iEt/h) / daK[HE,l]iua;uoO) 
[zTrh) Jo 

(93) 

with 

He = 4:u\Ho -E) = i^-Ae^- AEu ^ ) (94) 

\2m ) 

Notice that apart from an additive constant. He is just the harmonic oscil- 
lator hamiltonian in two dimensions. Thus we see that the HPIl in ( p3| ) is 
related to the propagator for the two dimensional oscillator and the relation 
is given by 

K\He, P = l]iua; uqO) = exp{i4:e^a/h)k°'^uia\uioO)k"'^U2cr\u2oO) (95) 

with k°'^'^ is the one dimensional oscillator propagator for mass fi = m and 
frequency uj^ = —8E/M. 

Step 4 ■ The HPI2 /C of (^) is thus known. We use analysis of Sec. 3 to 
relate HPI2 for Hq in the cartesian and the parabolic coordinates. For this 
we notice that x, y defined by (|8^) do not change if we substitute 



Ui —Ml, and U2 — > —U2- (96) 

Therefore, we have 

6{f-ro)=g-^'^[6{u-u^) + 6{u + u^)] (97) 

and 

/C[/7o,l,l](rt;n)0) = [/C[ifo,4u^4M=^](a;uoO) + /C[i^o,4u^4n^](n^;-^ioO)] 

(98) 
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Step 5: Thus the desired quantum mechanical propagator Kq is obtained 
from ( p2D and (^) using the expressions (^) and ( ]95|) . Writing the answer 
in terms of the energy dependent Green function we get 



Go(F,n)|E) = j exp(4zeV//i)[kf^(Mia|MioO)kf^(M2^|M2oO) 

+kr - MioO)kr(M2(T| - M2O)] (99) 

This result agrees with equation (28)- (30) of ref. 20 after the integration 
variable a is replaced by r = Aa. 

7. H-Atom in three dimensions 



7.1. K-S Transformation: The H-atom in three dimensions is related to 
isotropic harmonic oscillator in 4-dimensions via K-S transformation. Local 
time transformations in the path integral formalism and K-S transformation 
have been used to relate the Coulomb problem Green function to the Green 
function for the oscillator problem in 4-dimensions. This relation is derived 
in our canonical formalism of path integrations follows. 



Step 1: Consider the 4-dimensional oscillator problem defined by 

n 2 1 

™ = + -Mn^u ^ (100) 
^ 2M 2 ^ ^ 

where 



U = (ui,M2,M3, W4) (101) 
Pu = {Pu„Pu2,Pu„Pu^) (102) 

We introduce the new co-ordinates (xi, X2, xs, X4) , in two steps. In the first 
step the transformation equations 

. 9 .a + cl). ^ 6 . ,a — (p. 



Ui = v^sm-cos(^ — ^ — } ; ^2 = v cos - sm^ ^ ^ 
U3 = V'"cos-cos(— ^ — ) ; Ui = y/r sm - sm( ^ ) (103) 
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define the co-ordinates r, 6,(j),a . In the next step we identify X4 = a and the 
equations relating r, 9, to {xi, X2, xs) = x are taken to be the same as those 
giving the relation between the spherical polar to cartesian co-ordinates in 
three dimensions. Thus 



Xi = r sin 6' cos 0, X2 = r sin 6' sin 0, X3 = rcos6', X4 = a (104) 

The classical hamiltonian ( p.00| ) when expressed in terms of the co- ordinates 
{x,a) and conjugate momenta {p^Pa) becomes 

X1P2- X2Pi)zpa\ , 1 ,^^2 




-MVL^r (105) 



Also 



Mr{xl + xl) \ 2 



rf^M = 16r sin dr d9 dcj) da 

d'^u = [IQ / r)dxidx2dx2,da (106) 



V2 = (le/^) = J (107) 
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Step 2: Next we set up path integrals for the 4-dimensional oscillator 
problem in two sets of co-ordinates systems. In the cartesian co-ordinates u 
the relevant HPI2 is equal to the HPIl for the same hamiltonian 

K,[Hr, 1, 1] = K[Hr, 1] (108) 

The path integrals in (|108D at t = become equal to 5^{u — Uq). Thus two 
HPI's in ( |108| ) coincide with the exact quantum mechanical propagator, k^'^'^ 
for the four dimensional oscillator problem with mass M and frequency Vl. 
To quantize the oscillator ( |10U| ) in coordinates (x, a) we set up the HPI2 
1C[HI^^, J, J]{x, a, t; xq, ao, 0) which is given by 



/C[if4*'^, J, J] (x, a, t] xo, ao, 0) 

/d£ 
——— exp{—iEt/h) / daK[H£, p = l]{x, a, a; x,aoO) (109) 
[ZTTh) Jo 
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Taking inverse Fourier transform this relation is written as 



roo 

/ daKlHs, p = l\{x, a, a; X, aoO) 
Jo 

= — f dtexp(iSt/h))C[H2"'',J,J]{x,a,t-,xo,ao,0) (110) 
16 J 

The hamiltonian function appearing in the HPIl of r.h.s of (|109|) is given by 

He = J{Hr-S)/lQ (111) 
\2M Mr{xl + xl) 2 ' j ^ ^ 



A factor of 16 appears in (|109| ) to ( |111| ) because a trivial scaling of time has 
been included in ( |109|) 

Step 3: Next we note that two HPI2 set up for the oscillator problem, 
appearing in ( |108| ) and ( |109| ), are normalized differently. Their values at 
t = are related by 

{l/J)5\x - xo)S{a - ao) = [S\u - uo) + S\u + uq)] (113) 
Thus we must have 



/C[i7f ^ J, J] (f , a, t; xq, ao, 0) 
= /C[ifr, 1, l]{ut; UoO) + lC[Hr, 1, l\{ut- -UqO) 
= kr(Mt|MoO) + kr(Mt| - MoO) (114) 

where the oscillator propagator, 'k°^'^ appearing in the above equation is the 
four dimensional oscillator propagator. 

Step 4 : We next set up path integral quantization of Coulomb hamiltonian 
in cartesian co-ordinates. The required energy dependent Green function is 
given by (^) in arbitrary co-ordinates. In cartesian co-ordinates we have 

^oo 

Go(f,fo|^)=/ daK\HQ- E,p = l]ixa]XQQ) (115) 
Jo 
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where Hq is the Coulomb hamiltonian 

Ho = ^-- (116) 
2m r 

Step 5: We can now assemble the above results to obtain the H-atom 
Green function. It should be noticed that for = the function Hs is the 
Coulomb hamiltonian apart from an additive constant |Mfi^ .Therefore, the 
integral of K[H^ , 1] (x, at;xQ, ao 0) over Aq;(= a — oq) is related to the HPIl 
for the function Hq — E 



K[Ho- E,p=l]{xt;xoO) = d{Aa)K[H£,l]{x,at;xo,aoO) (117) 

Jo 

where we have identified m = M/4, S = and ^ = £ . The integral of 

the left hand side of ( |117| ) over t is just the energy dependent Green function 
( |115|) for the Coulomb problem . Therefore using ( |110|) and (|112|) we get, 



Go{x, Xo\E) = — / daexp{ie^a) / d{Aa)IC[H2^'', J, J]{x, aa; xq, aoO) 
16 J Jo 

(118) 



which on using (|114D gives 



Go{x,xo\E) 

]^ roc rAiTT 

= — daexp{ie'^a) c/(Aa)[A;f'^(Ma, -MqO) + kf '=(Ma, -MoO)](119) 
16 JO Jo 

This is the desired connection between the solutions for the four dimensional 



oscillator and the Coulomb problems. The above result ( |119| ) is in agreement 
with the equation (108) of ref 19. 

7.2 Parabolic Coordinates: In this subsection we shall show how exact 
path integration can be done for H-atom in parabolic coordinates within the 
canonical formalism. In this case we proceed directly to set up the path 
integral representation for the propagator and evaluate it. The coulomb 
hamiltonian 
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-» 2 2 

Hcoul = ^-- (120) 
2m r 



takes the form 



2m(C + ^) 2m(e + ??) 2m^r] ^ + 
in the parabohc coordinates ^,ri,(f) defined by 

X — ^r/cos0 
y = ^?7sin0 

^ = ^(^'-^') (122) 

Also 

dxdydz = ^r]{^'^ + rj'^)d^ dr] d(f) (123) 
We, therefore, set-up the HPIl /Cfif^, y/g] in the parabohc co-ordinates with 

He = g^'\H,oui-E) 

and use it to define HPI2 fC as follows 

^[H, ^/g^/g\{^r](t), t; Co?7o0o, 0) 

/dE 1"°^ 
-^eM-^Et/h) daK[HE,mv^,t;^oVo(l>0,0) (125) 

/C satisfies the Schrodinger equation appropriate to the Coulomb potential 
but has the initial value 

\i^lC[H, ^/g,^/g\{iv <Pt; Co 0o 0) = ^^^-^^^^:^<^(C-6)5(^-^o)5(0-0o) 

(126) 

whereas the Coulomb propagator should be normalized to 
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-1 oo 

Si^-^o) = -^7^—^J{^-^o)S{v-Vo) E '5(0-00 + 27™) (127) 

><V I'b ' V ) n=—oo 

Taking ( |127[ ) it into account the Coulomb propagator will be given by 



Ko(ft|roO) 

oo 

= E /C[ifo,v^,v^](e# + 27rnt;eor/o0o,O) (128) 

n=—oo 

dE 



——exp{-iEt/h) / da ^ ^[^^i?, 1](^'7, + 27rn, a; ^0^7000, 0) 

(129) 



Using time integrated version of (pTl) and expressing change in scaling as 
change in potential for HPl we get, with a = C,ti , 

CO roo 

daK[HE, l](e^0cT; ^oVoM = / K[hE, l](e^0a; CoVo4>oO)da (130) 

JO 

where 



p| / l 1 \ 2 2 2 ^ ( ^ ^\ 

^ 2m ~^ 2m ~^ 2m ~^ ) ^ ^^2 ~^ ^2 J ^ ) 

Thus we have 

dE 



^—^^e^^{-tEt/h)Go{r,fo\E) (132) 



with 



/■oo °° 

Go{r,ro\E)= da K[hE,m,vA + '27in,a-,Co,Vo,(po,0) (133) 

•^0 n=-oo 

We shall now write the canonical path integral K[hE,l] in the discrete from 
and do the and path integration explicitly. For this purpose we write 
K[hE, 1] of (jisg) as 



28 



„iV-l N-l \ 

/ n (dijdrjj) Yl dp^jdp] exp 

J j=l \j=o / 



. k k 

(134) 



exp 



8m 



1^1 
Cfe Vk. 



(135) 



where 



N-l ( 
3=0 \ 



'j Wj+i 



Pi 

2m VC| ' Vj. 



(136) 

The (pk integral may first be done giving rise to the product of delta functions 
5{pf. — Pk+i) which may in turn be used to do all except one pf. integration. 
The last remaining p'^ integral is Gaussian and is easily written down. This 
gives 



= {Ih/2mT) exp ( ^ 



where 



Using the identity 



t (—) 



1/2 



I exp 



1 oo 



'il{(f)- (f)o + 2Trny 
2Th 



(137) 



(138) 



oo / X \ 
n=-oo ^'^^^ 

We get the Coulomb Green function as 



(139) 
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Go(r-',fo|E)= / i^[^^„,l](e,^,t^;eo,r/o,O)expH0-0o)]c?or (140) 



where Hn in ([L40|) is the hamiltonian function 



2m 2m ' ^ 2m^ a'^^^ t]^' ^ ' 

= H{^,p^) + H{r],p,) (142) 



As Hn is a sum of two hamiltonians each depending only on one set of 
conjugate variables, we have 

K[Hr„ 1] (e V cr; Vo 0) = K[H, 1] (e^o; ^0)^^ 1] (r/ao; r^O) (143) 
The HPIl for H is known from Sec. 5.1 and hence we get 

K[HnA]=(^^^^](^^^^]jJ (144) 
\^ zh sm ut J yzn sm cut y \zh sm tut y \^ sm ut J 

exp I ^ cot ut{e + eo + Vo + Vo) ] (145) 

oo 

K = — exp[tni<p-<Po)]K[Hn,m,v,t;^o,Vo,0) (146) 
The required Coulomb Green function is then given by 

= / da ey.-pi—2ie^a ITi) (— ) exp(m(0 — 0o)) 

JO ^rr'oo V^'^sincuo-y 

f ^"^'^ r^9 ^2 2 2nl T 

exp|^^cot^(T[^- + ^o +r/ +r]Q\\Jn 



" sin Lua 



jJ'^™] (147) 

V nsmuja 
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Using the identity 



J2 exp{ina)Jn{r)Jn{p) = MR) (148) 

n=—oo 

where 

= r^ ^p"^ + 2rpcosa (149) 



the Green function (|147|) becomes 



Go(e,^,0;eo,r?o,0o|^) (150) 

= r da exp(2.e V/;.) (-^) ' Jo f , r] 
Jo \insmua / \nsmua ) 

exp|^cota;t(e2 + eo + r?' + %')} (151) 



with R given by 



i?^ = + + 2^^o cos(0 - 0o) 

= rroCos(V'/2) (152) 

and 

cosV" = cos 6* cos 6^0 + sin 6* sin 6^0 cos(0 — 0o)- (153) 

Defining k = {muj / ih) the above result is seen to coincide with that obtained 
by Ho and Inomata in ref 13. 

7. 3 Radial propagator for Hydrogen atom : The hamiltonian for the radial 
hydrogen atom problem is 

// = ^-£! + 5!i(i±i) (154) 

Note that 1/r^ term in the hamiltonian can be thought of as coming from 
a pg/2mr^ term in hamiltonian for corresponding two dimensional problem. 
Therefore, this path integral solution of this problem can be obtained by 
relating it to the Coulomb problem in two dimensions. This relation can be 
derived in a way parallel to the derivation in Sec. 5. 
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8. Concluding Remarks 



We have summarized the main properties of the hamiltonian path integral 
formalism of rcf 16, 17, and 18. We have also discussed application to exact 
path integration. Here, we wish to highlight the following points. 

(a) For all other existing literature on path integration, local scaling of 
time is a technique used in exact solution, where as for us it is an essential 
ingredient of quantization scheme itself. In some cases, such as in the example 
of H- atom in Sec. 6 , no further local scaling of time is necessary to complete 
exact solution. 

(b) The best feature of this approach is that path integral quantization 
is formulated in terms of the classical hamiltonian in arbitrary co-ordinates. 
Performing changes of variables becomes quite trivial in the path integral 
formalism in our scheme, as one can do this at the classical level itself and use 
the classical hamiltonian in arbitrary coordinates to set up the canonical path 
integral quantization. The boundary conditions are correctly reproduced 
when initial condition for the propagator is correctly taken into account as 
discussed in Sec 3. In this respect our approach is new and more direct than 
all other existing treatments. 

(c) The canonical approach offers more general possibilities for local scal- 
ing of time, than those exploited here. It should be possible to generalize the 
scaling relation to the case when f{q) is replaced by a function of p or by a 
function of p and q. It is not known whether this can be used to do exact 
path integration of new problems. 

(d) A class of potentials, such as Scarf potential, Rosen Morse potential 
etc., can be related to a free particle on sphere. Though we are not able 
to do exact path integration for free particle on a sphere, our scheme can 
be used to obtain inter-relations of propagators for the potentials within the 
class. A possible approach to solve particle on sphere problem exactly may 
be to formulate it as a path integral for a constrained systems. For this a 
discretized form of path integral for the constrained systems is needed. 

(e) In ref. 18 correspondence between the canonical and lagrangian forms 
of path integrals was estabhshed. For most of the problems, not taken up for 
discussion in this paper one can, in principle, utilize correspondence between 
the hamiltonian and the lagrangian forms of path integration and follow the 
existing literature to complete the path integration. However, this is not very 
interesting. 
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(f) Doing a close analysis of discussion in sections 5,6 and 7, we see that 
with exception of the path integral in Sec. 7. 2, hardly any path integral, is 
actually computed. All the details of computing and manipulating path 
integrals are buried in the results summarized in Sec. 2 which were derived 
in ref. 18. The present paper thus isolates and identifies a set of properties 
of Hamiltonian path integrals useful for path integration. 

(g) The aim of path integration formalism is not merely to solve those 
quantum mechanical problems exactly for which the exact solution can al- 
ready be obtained by several methods. This activity is useful only if it gen- 
erates new techniques or gives new insights into the path integral formalism. 
It would be far more interesting and useful to explore the generalizations and 
advantages of the hamiltonian approach for path integration. 

(h) For the constrained hamiltonian systems, in a most direct but a for- 
mal approach, one sets up path integral over the reduced phase in terms of an 
independent (reduced) set of conjugate variables {q*,p*) which, in general, 
are functions of an initial set {q,p). By means of a formal canonical transfor- 
mation, when the path integral is expressed as path integral over the initial 
set {q,p) over full phase space, Dirac delta functions enforcing constraints 
and determinants appear in the integrand. This path integral requires a pre- 
cise lattice definition. It then appears that, any attempt to precisely define 
a discrete version of path integration for the constrained systems must solve 
the question of incorporating canonical transformation as change of variables 
in the phase space path integral approach. This merits further investigation. 
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